1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 300,000 iteration;
- 200,000 burn-in;
- 100 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 0.999 1.00
## LDevsum 0.999 1.00
## dh0 1.015 1.05
## dh1 1.015 1.05
## dl0 1.868 3.52
## dl1 1.718 3.22
## dl2 1.336 2.05
## dl3 1.096 1.23
##
## Multivariate psrf
##
## 1.78
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1224.97003 | 22727.5825 |
| DIC3 | 1161.88355 | 20648.2725 |
| PWAIC | 44.02095 | 996.1293 |
| WAIC | 1189.68711 | 21027.5778 |
2. Separate Modeling
2.1 Trace plots for convergence check
The current MCMC setting is:
- 300,000 iteration;
- 200,000 burn-in;
- 100 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.01 1.03
## LDevsum 1.01 1.03
## dh0 1.60 2.91
## dh1 1.52 2.73
## dl0 1.37 2.49
## dl1 2.57 5.95
## dl2 1.02 1.05
## dl3 2.43 6.74
##
## Multivariate psrf
##
## 1.87
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1448.01985 | 23921.990 |
| DIC3 | 1292.89309 | 21168.527 |
| PWAIC | 94.59621 | 1536.806 |
| WAIC | 1360.96116 | 21915.208 |